yunif <- array(); a <- 5;b <- 15
for(i in 1:10) yunif[i] <- ran.unif(a,b)
hist(yunif,prob=TRUE); title(sub="Uniform in 5,15")
yexp <- array(); rate <- 1/10
for(i in 1:10) yexp[i] <- ran.exp(rate)
hist(yexp,prob=TRUE); title(sub="Exponential of rate 0.1")
ynorm <- array();
for(i in 1:10) ynorm[i] <- ran.norm()
hist(ynorm,prob=TRUE); title(sub="Std Normal")
Run the code above in your browser using DataLab